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Abstract 



Aims. We provide simple and accurate analytic approximations for the low frequency double Compton emission coefficient that are 
applicable in a broad range of physical situations up to mildly relativistic temperatures. These approximations may be useful for 
checking in which circumstances the double Compton process is important. 

Methods. We perform series expansions of the double Compton emission integral for low energies of the incident photon and electron 
and compare the derived analytic expressions with the results obtained by numerical integrations of the full double Compton cross 
section. 

Results. We explicitly derived analytic approximations to the low frequency double Compton emission coefficient for initial 
monochromatic photons and Wien spectra. We show that combining the analytic approximations given in this paper, an accuracy 
of better than 5 % over a very broad range of temperatures and under various physical conditions can be achieved. The double 
Compton emissivity strongly depends on the ratio of the incoming photon and electron energies: for hard photons and cold electrons 
the emission is strongly suppressed compared to the case of similar photon and electron energy, whereas in the opposite situation, i.e. 
hot electrons and soft initial photons, the emission is enhanced. For photons and electrons close to thermodynamic equilibrium the 
double Compton emissivity increases less rapidly with temperature than in the Lightman-Thorne approximation and the corrections 
exceed ~ 10% for temperatures above 4keV. 

Key words. Radiation mechanisms: double Compton scattering — thermal plasmas 



1. Introduction 

Double Compton (DC) scattering is a process of the form 
e + yo «-» e' + j\ + j2- It is related to Compton scattering, 
like Bremsstrahlung is related to Coulomb scattering of elec- 
trons by heavy ions. It corresponds to the lowest order cor- 
rection in the fine structure constant a to Compton scattering, 
with one additional photon in the outgoing channel, and like 
thermal Bremsstrahlung it exhibits an infrared divergence at 
low frequencies. In spite of these similarities, in the case of 
Bremsstrahlung the cross section directly depends on the veloc- 
ity of the electron relative to the charged particle: there is no 
Bremsstrahlung emission for resting electrons (and ions). On the 
other hand, for DC scattering the dependence on the velocity of 
the electron enters only indirectly and can be obtained as a result 
of special relativistic coordinate transformations. Therefore DC 
emission occurs whenever there are photons and free electrons, 
no matter what the temperature of the electron gas is. 



It was realized by iLightmanl (Il98ll). iThomd (1 198ll) and 
later by iPozdnyakov et al.l d 19831) and ISvenssonl (11984 that the 



DC process in comparison to Bremsstrahlung can become the 
main source of soft photons in astrophysical plasmas with low 
baryon density and in which magnetic fields are negligible. It 
has been shown that for given number densities of protons, 
N v , and photons, N y , the DC process is expected to dominate 



over Bremsstrahlung wherQ/Yp < lON y (kTJm e c 2 ) 5/2 is ful- 
filled. This implies that due to the large entropy of the Universe 
(there are ~ 1.6 x 10 9 photons per baryon), at early stages (red- 
shifts z > 5 x 10 5 , i.e. at temperatures higher than ~ 10 6 K), 
the DC process becomes very important for the thermalization 
of possible spectral distortions of the cosmic microwave back- 
ground (CMB) and the evolution of chemical potential distor- 
tions after any significant release of energy in the early Univer se 
(ISunv aev & Zeldovich, 1970; Illarionov & Sunvaevl Il975alfbt 
Danese & de Zotti, Il982t iBurigana. Danese & de Zottil Il99lt 
iHu & SilkLll993l) . Other possible environments in which the DC 
process could be of some importance might be found inside the 
sources of y-ray bursts or two temperature accretion disks in the 
vicinity of black holes, both in close binary systems and active 
galactic nuclei. Given the relevance of DC emission as a source 
of soft photons within the context of the thermalization of CMB 
spectral distortions and hard X-ray sources, it is important to in- 
vestigate the validity of the approximations usually applied to 
describe the DC process. 

Expressions describing DC emission in an isotropic, cold 
plasma and for low energy initial photons, i.e. small loq = 
hvo/m e c 2 , where vo i s the frequency o f the initial photon yo, 
were first deri ved by ILightmanl dl98ll) and independently by 
iThornd dl98ll) within the soft photon limit. In this approxi- 
mation it is assumed that one of the outgoing photons {j\ or 
72) is soft compared to the other. Under these assumptions 
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1 For this estimate it was assumed that the incident photons are 
distributed according to a Wien spectrum with effective temperature 
T c = T y and that the plasma is neutral, N e = N p . 
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the DC cross section in creases oc aj 2 with photon energy (see 
iJauch & Rohrlichll 19761 Eq. 11-45) and the DC emissivity for 
broad initial photon distributions is proportional to the integral 
/ ~ J Vg«(vo)[l + n(vp)] dvn , wher e here n(v) is the photon oc- 
cupation number (Lightman, 1981). We give more details about 
the Lightman-Thorne approximation in Sect. 12.21 

Again presuming cold electrons and low energy incident 
photons but ma king a priori no assumption about the emitted 
photon energy, iGouldl (Q984) obtained an analytic correction 
factor for the DC emission coefficient of monoenergetic initial 
photons, which extends the Lightman-Thorne approximation be- 
yond the soft photon limit. Here the corrections are directly re- 
lated to the increase of the emitted photon energy (up to energies 
~ Vn/2), while the low frequency emission spectrum remains un- 
changed. 

However, when the temperature of the electrons increases 
or when the energy of the incident photons grows, one expects 
corrections to become important, an aspect that has not been in- 
cluded in either of the aforementioned approaches. The correc- 
tions due to motion of the electron can be obtained starting with 
a resting electron in an anisotropic radiation field and then per- 
forming the corresponding transformations into the frame where 
the electron is moving. In contrast to this, the dependence on the 
energy of the incident photon is related to the exact formulae for 
the DC cross section as computed using the general theory of 
quantum-electrodynamics. Here we treat the electron-photon in- 
teraction due to DC scattering quantum-electrodynamically for 
moving electrons and perform all the computations directly in 
the lab frame, where the photon and electron distributions are 
assumed to be isotropic. 

Temperature corrections to the low frequency DC emissivity 
for an incoming Wien sp ectrum, again wi thin the soft photon 
limit, were discussed by [Svenssoa {i984). It was shown that 
in the mildly relativistic case the low frequency DC emissiv- 
ity increases significantly slower with tempe rature than in the 
Lightman-Thorne approximation. However, ISvenssonl (l984) 
only treated one specific case and a more general extension of 
the Lightman-Thorne approximation is still missing. 

In this paper we wish to extend the Lightman-Thorne ap- 
proach to cases hvo < m e c 2 and kT e < m e c 2 for isotropic ini- 
tial photon distributions, but still within the soft photon limit. In 
addition, we will focus on situations when DC absorption and 
stimulated emission are negligible. The extension to cases when 
(l>2 ~ coo and a discussion of the full kinetic equation for the DC 
process will be left for a future work. We study the DC emis- 
sion integral both numerically and analytically and derive ap- 
proximations to the DC emission coefficient, which are valid in 
a very broad range of physical situations up to mildly relativistic 
temperatures. We show that the DC emission coefficient can be 
expressed in terms of the Lightman-Thorne emission coefficient 
times a corrections factor, which can be regarded as very simi- 
lar to the Bremsstrahlung Gaunt-factor, although it has a com- 
pletely different physical nature. As an example, in Fig. [4] we 
summarize the results for this low frequency DC correction fac- 
tor in the case of monochromatic initial photons. One can see 
that the Lightman-Thorne approximation is accurate in a very 
limited range of photon energies and electron temperatures. We 
also investigate the range of applicability of our approximations 
for initial Wien spectra (Sect. 1431) . 

In the following we use h-c-k- 1 . 



2. DC emission for monochromatic photons and 
thermal electrons in the soft photon limit 

In this Section we give a detailed derivation of the approxima- 
tions describing the DC emission for monochromatic initial pho- 
tons and thermal electrons within the soft photon limit. In par- 
ticular, we provide some additional comments on the Lightman- 
Thorne approximation (Sect. l2~2l i and introduce a DC correction 
factor (Sect. l2.3l l relative to the Lightman-Thorne formula which 
will be used in the following. See Table[T]for overview. 

2.1. General formulation 

DC emission is a result of the interaction between an electron 
and a photon with one additional photon in the outgoing channel: 

e(P) + y(K ) — » e(P') + y(K x ) + y(K 2 ) . 

Here P = (E,p), P' = (E',p') and K t = (v,,Jt,) denote the 
corresponding particle four-momentcQ. The full DC scattering 
squared matrix element and differenti al cross sections for vari- 
ous limiting cases were fi rst derived bvlMandl & Skvrmd(ll952l) 
and may also be found in lJauch & Rohrlichl(ll976l) . 

In this Section we focus our analysis only on the emission 
process for isotropic, monochromatic initial photons, with phase 
space density n(v) = A^o S(v— Vo)/vSL and isotropic, thermal elec- 
trons of temperature T e . Above No in physical units is defined by 
N y $ = SnNo/c 3 , where N y fl is the photon number density. In this 
case the change of the photon phase space density, n 2 = n(v 2 ), at 
frequency v 2 due to DC emission can be written as 

dn 2 m No r 3 r r do- DC 

ir^rTfJ dp J dQo J (1) 

Here the DC differential cross section (cf. IJauch & RohrlichL 
119761 Eq. 11-38) is given by 

do-pc _ u r l vi vi X 

dA ~ (An) 2 g^yvo (P + K - K 2 ) ■ K x ' 

where X may be found in the Appendix [B] and we defined the 
Mpller relative speed of the incident electron and photon as 
g0i — P ■ K, ■ — 1 -f3p e j, with the dimensionless electron velocity 
/? = |v|/c = \p\ - \p\/E and the directional cosing p e , - p ■ k,. 
Furthermore a = e 2 /4n « 1/137 is the fine structure constant, 
ro = a/m e « 2.82 x 1CT 13 cm is the classical electron radius and 
y = y/l - p 2 denotes the Lorentz-factor of the initial electron. 

In Eq. (Q]) we have neglected induced effects, i.e. we assumed 
No <k Vy and that the electrons are non-degenerate. The lat- 
ter simplification, for sufficiently low temperatures and electron 
densities (such as in the Universe for T e < m e ), is justified, but 
as we will discuss in Sect. 13.21 for broad initial photon distribu- 
tions, stimulated DC emission can become important. However, 
a full treatment of this problem is much more complicated and is 
beyond the scope of this paper. The electron phase space density 
may be described by a relativistic Maxwell-Boltzmann distribu- 
tion 

f(E) = — e- EM , (3) 

4nmlK 2 (l/e e )0 e 

2 Henceforth bold font denotes 3-vectors. 

3 Note that in the following an additional hat above 3- and 4- vec- 
tors indicates that they are normalized to the time-like coordinate of the 
corresponding 4- vector. 
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wher e A~2(l/0 e ) is the modified B essel function of the second 
kind (iAbramovitz & Stegun|jl965]) , with 8 e = T e /m e , and where 
N e is the electron number density, such that = J f(E)d 3 p. 
In the low temperature limit (6 e <<c 1) the relativistic Maxwell 
Boltzmann distribution © can be handled by the expression 
given in Appendix [A] which is useful for both analytic and nu- 
merical purposes. 

In general the DC emission integral (Q} has to be performed 
numerically. Some comments on the numerical approach we use 
to solve the multidimensional Boltzmann integrals is given in 
Appendix ICl However, in the limit of low temperatures and en- 
ergies of the initial photon it is possible to derive various use- 
ful analytic approximations, which we discuss in the follow- 
ing. In practice we obtained the numerical results for the soft 
photon limit by setting the frequency of the emitted photon to 
a very small value as compared to the incident photon energy 
(v 2 ~ v x 10~ 4 ). 



2.2. Lightman-Thorne approximation 

The expression for DC emission fro m monoenergetic initial pho- 
tons and col d electrons as g iven by iLightmanl (1 1 98 lb and inde- 
pendently bv lThornel (1 19811) can be deduced from the emission 
integral (HJ by performing a series expansion of the DC differen- 
tial cross section © in lowest orders of vo «: m e and v 2 <s m e 
and setting ft — 0, i.e. assuming that the electrons are initially at 
rest. In this limit f(p) = N e 5{p)IAnp 2 and the integration over 
d 3 p can be carried out immediately. 

The series expansion in lowest order of v 2 <s m e is equiv- 
alent to using the soft photon approximation for the DC differ- 
ential cross section, i.e. assuming in addition that V2 <s vo. Due 
to energetic arguments the scattered photon frequency has to be 
close to the initial photon frequency (vi ~ vo). Therefore the 
infrared divergence due to the integration over the phase space 
volume for the photon y\ is automatically avoided. In the follow- 
ing we use the notations pij = £, • kj and <p,j for the directional 
cosines and azimuthal angles between the photons i and j, re- 
spectively. It then follows: 



j soft 
uu DC 



dA 



/ ,2 

Ct (1 1 



1 - 



A/j 2 



(4) 



with A// = [iQ2 - fin- Furthermore we have introduced the di- 
mensionless photon energy <y, = v,/m e . 

Now, after aligning the z-axis with the direction of the emit- 
ted photon all the integrations can be easily performed analyti- 
cally. Making use of the identity yuoi = ^O2Mi2+cos(0o2— 0i2)(l- 
^02 ) (1 ~/Ai) l ^ 2 > wnere 002 and <fiyi denote the azimuthal angles 
of the initial and the scattered photon, respectively, this leads to 
the Lightman-Thorne approximation for the DC emission spec- 
trum of cold electrons and soft, monochromatic initial photons 



<9«2 



em,L 



4a 
3^ 



N e No <X T ■ 



(5) 



with the Thomson scattering cross section cr T = 87rr ( ,/3 « 
6.65 x l(T 25 cm 2 . The expression (0 is the DC scattering equiv- 
alent of the Kramers formula for thermal Bremsstrahlung. The 
characteristics of equation 01 can be summarized as follows: 
the DC emission spectrum increases <x a> 2 and, as mentioned 
earlier, exhibits an infrared divergence for V2 — > 0. It is usually 
assumed that under physical conditions inside an astrophysical 
plasma, emission and absorption balance each other below some 



Table 1. Analytic approximations for the DC correction fac- 
tor (0 relative to the Lightman-Thorne formula as obtained in 
Sect. |2] Within the quoted range of applicability each approx- 
imation should be accurate to better than ~ 5%. However, for 
G™(w , 6» e ) = (G'7) th as based on Eq. (QUi and tO we refer 
the reader to Fig. 



Symbol 

erVo) 

G»(0e) 

G£(0 e ) 
G m (QJo,6 e ) 



Reference 
Eq. ©,Fig.[T] 
Eq. ©, Fig.0 
Eq. (Ell, Fig. E] 

Eq. tTM . Fig.[2] 
Eq. (B}, Fig. [3] 

see above, Fig. [3] 



Range 



comment 



6, =0,a> < 0.15 
£) c = 0, oj < 1 
9 e < l,w « 1 
6) c < 0.5, w a 1 
G e < 0.2, co < 0.1 
Fig.Sl 



expansion 
inv. factor 
expansion 
expansion 
expansion 
inv. factor 



frequency, V2, m in, and that energy conservation in addition intro- 
duces some high frequency cutoff, V2, ma x- Therefore the photon 
production rate increases logarithmically, 



dN 2 
~dt 



em,L 



4ff 2 , / V 2,max 

— N e N cr T co z x In 

in \ V 2 min 



(6) 



with the ratio of the upper to the lower frequency cutoff of the 
DC emission spectrum. 

It is obvious that due to the DC process the interacting elec- 
trons and high frequency photons can lose some part of their 
energy due to the emission of soft photons. However, in the ma- 
jority of applications the energy losses connected with Compton 
heating and cooling are much larger. 

As will be discussed in more detail below, the Lightman- 
Thorne approximation only describes the low frequency part of 
the DC emission spectrum for cold electrons and soft initial pho- 
tons (vo/m e < 10~ 3 ) well. For initial photons with slightly higher 
energies, next order terms lead to a significant suppression of the 
emission at low frequencies relative to the Lightman-Thorne ap- 
proximation (see Fig. [T]and the discussion in Sect. 12. A\ . 



2.3. DC correction factor within the soft photon limit 

Below we discuss the numerical and analytical results for the 
DC emission integral ([]]). For convenience using Eq. (0 and the 
Lightman-Thorne approximation we introduce the new func- 
tion G as 



G m (a>o,e e ) = 



d,n 2 \l 
d,n 2 



(7) 



em,L 



in order to compare the different cases. This function can be re- 
garded as very similar to the Bremsstrahlung Gaunt factor, al- 
though it has a different physical origin. Deviations of G from 
unity are directly related to DC higher order corrections in the 
energies of the initial photons and electrons. 

2.4. Cold electrons and energetic initial photons 

We expand the DC differential cross section (fSJ in the lowest or- 
der of u>2 = v 2 /m & and setting fj = 0. Then we expand the result- 
ing expression up to 4th order in ojq to take into account higher 
order corrections in the energy of the initial photon. Carrying out 
all the integrations in the emission integral yields the correc- 
tion factor G^((Wo) following from the direct series expansion: 



G m (co ) = 1- —oj + — co Q - 



357 2 7618 3 21498 



175 



0> n H Wn . 

175 



(8) 
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Figure 1. DC emission correction factor G Q m for cold electrons 
and monochromatic initial photons as a function of cjq = vo/m e . 
Also shown are the direct expansion ([8]), up to different orders 
in wo, and the inverse approximation (0. For comparison, we 
also show the frequency dependence of the total Compton cross 
section, ctkn, as given by the Klein-Nishina formula and the cor- 
rection factor k(u)q) = G[J 1 /(cr KN /cr T ) per act of Compton scat- 
tering. 



Figure [TJ shows the full numerical result for G^cjq) in compar- 
ison to the analytic approximation (0, taking into account the 
corrections up to different orders in coq- The approximation con- 
verges only very slowly and in the highest order considered here 
it breaks down close to a>o ~ 0.15. Due to this behavior of 
the asymptotic expansion one expects no significant improve- 
ment when going to higher orders in a>o, but the monotonic de- 
crease of the emission coefficient suggests that a functional form 
GjJ, = [1 + Yj\=\ a k(Jo\ l could lead to a better performance. 
Determining the coefficients a,- by comparison with the direct 
expansion (O one may obtain an inverse approximation for the 
DC correction factor 



GS; m Vo) : 



l 



1 + 



o» 



84 2 
25 % 



2041 

875 



9663 4 ' 
4375 % 



(9) 



As Fig. [TJ clearly shows, Gm inv («o) provides an excellent de- 
scription of the numerical result up to very high energies of 
the initial photon. For comparison, we also show the frequency 
dependence of the total Compton cross section, cr KN , as given 
by the Klein-Nishina formula and the correction factor k(u>o) = 
G^/(cr KN /cr T ) per act of Compton scattering. Since it is ex- 
pected that per single Compton scattering a small fraction of 
electrons and photons undergo DC scattering, this comparison 
hints towards the fact that a large part (~ 2/3) of the decrease in 
the DC emissivity with u>q is due to the decrease in the proba- 
bility of Compton scattering photons with larger energies due to 
quantum-electrody namic corrections . 

2.5. Hot electrons and soft initial photons 

We now consider the limit vo <K m e and v 2 <K vo. Expanding the 
DC differential cross section in the lowest orders of vq and V2 
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Figure 2. DC emission correction factor G™ for soft initial pho- 
tons (wo = 1CT 4 ) as a function of the electron temperature 
0e = T e /m e . Also shown is the full integral approximation ( 112al ) 
and the expansion ( I12bl ). taking into account the corrections up 
to different orders in 6* e , as indicated. 



one obtains 

j soft 
uu DC 



dA 
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2^2 



2 2 

g0O v 2 



-y 2 A A l (y 2 A A l 



y 2 /lo/li^a;oi 



-a i)4 
^oari2) 2 



'01 



(10) 



with the electron and photon scalars Aj = P ■ K \ = 1-/3 yU e / 
and a,j — K, ■ Kj — 1 — /j/j. Now we perform all the angular 
integrations of the Boltzmann emission integral ([TJ but assume 
monoenergetic electrons, with velocity /3o- This then leads to 



dn 2 



<9«2 

~dt 



(11) 



em,L 



Here the factor of y? can be immediately understood since the 
initial photons inside the rest frame of the electron will on aver- 
age have an energy of af^ ~ jqgjq and hence the DC emissivity, 
i.e. oc (oJq) 2 , will be ~ j\ times larger than for resting electrons. 

From Eq. ( fTTT i the correction factor for thermal electrons 
with temperature 6> e can be found by averaging over the rela- 
tivistic Maxwell-Boltzmann distribution (f3j): 



r 



f(E )p 2 dpo=2{y 1 ) 



1 



[i +24e 2 ]K (i/e e ) + 8e e [i + 6(? 2 ] Ki(\/o e ) 



K 2 (l/O e ) 



(12a) 



- P 2 y l/2 . Here 



with Eq = y m e , p = y m e /3 and y = (1 - Ho 
Ki(l/6 e ) denotes a modified Bessel function of kind i. For low 
temperatures one can find the approximation 



G™(6 e ) « 1 +66> e + 15f^ + 



45 



(12b) 



45 

Figure [2] shows the numerical result for G™ in comparison with 
the analytic approximations (fT2l . The low frequency DC emis- 
sion strongly increases with temperature. One can gain a factor 
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Figure 3. DC emission correction factor G m for different ener- 
gies of the initial photons as a function of 8 e = T e /m e . Also 
shown is the full integral approximation (fT4l in combination 
with the inverse factor (I16ab and, for some cases, the expansion 
(l5[ taking into account correction terms up to fourth order, as 
indicated. 



of a few for 6> e < 0.5 and even a factor ~ 30 for 8 e ~ 1. The 
expansion ( I12bl ) provides an excellent description of the numer- 
ical results up to high temperatures {8 e ~ 0.5), but for higher 
temperatures the integral approximation d!2ab should be used. 

2.6. Hot electrons and energetic initial photons 

To improve the analytic description of the DC emission we again 
return to the soft photon limit, expanding the DC differential 
cross section (J3J in the lowest orders of V2 <k m e . To include 
higher order corrections we also expand this expression up to 
4th order in cjq. We then carry out all the angular integrations 
of the Boltzmann emission integral ([T]) for monoenergetic elec- 
trons, with velocity fio. This leads to 



G m (« ,A)) = To 



l+$-yfl+24 + i$) ro wo 



357 
'^25 
7618 



(13) 



As in the previous case the correction factor for thermal elec- 
trons with temperature 8 e can be found by performing the 1- 
dimensional integral 

(G m \ h = G m (co o ,0 e )= G m (coo,/3 )f(E )p 2 dp . (14) 
Jo 

Unfortunately here no full solution in terms of simple elemen- 
tary functions can be given, but the integral can be easily eval- 
uated numerically. In the limit of low temperatures one can find 



the simple approximation 
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(15) 



for (O, which clearly shows the connection to the two limits 
discussed above, i.e. Eqn. © and (112bl >. 

As was shown in Sect. 12.41 an inverse Ansatz for G m led to 
an excellent approximation for the full numerical result in the 
cold plasma case. To improve the performance of the obtained 
analytic approximation for G m (o>o,j8o) we tr i e( l man Y different 
functional forms, always comparing with the direct expansion 
(foi l. After many attempts we found that 



G2 v (co ,/3 ) = 



1 + 2X1 /*W>) rK 



with the functions fkifio) 



ZiOSb) = 

fz(fio) = 
/3O80) = 
/4O80) = 



1 



1+^ 
1 
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(16b) 
(16c) 
(16d) 
(16e) 



provides the best description of the full numerical results. 

Figure [3] shows the numerical results for G m (u>o, 8 e ) in com- 
parison with the integral approximation (fT4l in combination with 
the inverse factor ( I16ab and the direct expansion (15[ taking into 
account correction terms up to fourth order. The performance of 
the integral approximation is excellent in the full range of con- 
sidered cases, but the direct expansion breaks down at lower and 
lower temperatures once the initial photon frequency increases. 
For wo < 0.1 the direct expansion should be applicable within a 
few percent accuracy up to 6* e ~ 0.2. 

In Fig|4]the dependence of the DC emission correction fac- 
tor on the energy of the initial photons and the temperature 
of the electrons is illustrated. From this one can see that the 
Lightman-Thorne approximation is applicable in a very limited 
range of photon energies and electron temperatures. On the other 
hand the approximation based on ( fToT l works in practically the 
full considered range with very high accuracy. Focusing on the 
curve G m = 1 one can see that for low energies of the initial 
photon (0)0 < 10~ 2 ) the required electron temperature scales 
like 9 e ~ 0.7 u>o- Considering the first order correction terms 
in Eq. (fT3T l one also finds this scaling. Up to coq > 0.15 the 
necessary temperature of the electrons is 8 S < uq. This shows 
that moving electrons can easily compensate the suppression of 
DC emission due to the increasing energy of the initial photon. 
However, for wo > 0.15 the electrons have to be up to several 
times hotter to accomplish this. 
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0.511 




Figure 4. DC emission correction factor for monoenergetic ini- 
tial photons and thermal electrons. The lines indicate the con- 
tours of G m = const as obtained by full numerical integrations. 
Within the shaded regions the approximation based on ( TToT i is 
accurate to better than 1% and 5%, respectively. 



3.1. The effective DC emission correction factor 

Here we will give analytic expressions for the effective DC emis- 
sion correction factor with three different approaches, each with 
their advantages and disadvantages. We will discuss the range of 
applicability for these expressions for Wien spectra in the next 
Section. 



3.1 .1 . Approximation based on G™ 

As was shown in Section l2~6l the expression ( TToT l averaged over 
a thermal electron distribution (compare with Eq. (14[ ) provides 
an excellent description of the DC emissivity for monochromatic 
initial photons. If induced effects are negligible one can simply 
convolve this result with the corresponding initial photon distri- 
bution n( v) to obtain: 



dx . 



(20) 



Here ^G™ v ^ h denotes the thermally averaged expression ( TToT ). As 
will be shown below, the approximation d20b works very well in 
a broad range of different temperatures. However, it involves a 
2-dimensional integral which in numerical applications may be 
too demanding. 



3. DC emission for more general initial photon 
distributions 

Based on the results given in the previous Section one can obtain 
simple expressions for more general initial photon distributions. 
For initial photons with phase space distribution n(v) within the 
Lightman-Thorne approximation one finds 



dn 2 
~dt~ 



soft 



DC,. 



4a fit 

= N ? °t ~k x 7 o , 
3n xl 



(17) 



where here we defined the dimensionless photon frequency jc, = 
cJi/Oy = Vi/T T We also introduced the effective temperature T 7 
which characterizes the total energy density of the photon field. 
For Planckian photons T y is identical with the thermodynamic 
temperature. 

The DC emi ssion factor Iq is given by the integral (compare 
|Lightmaii<11981|>,Eq. (lOb)forn(v) « 1): 



/o 



Jo 



x n(x) dx . 



(18) 



over the initial photon distribution n(v). Stimulated DC emis- 
sion was neglected, but we will discuss possible extensions in 
Sect. 13.21 Here one factor of x 2 arises from the conversion 
iVo — > v 2 n(v), and the other is due to the scaling of the Lightman- 
Thorne approximation Eq. (0) with coo. Note that here and below 
we omit the index '0' for the frequency of the initial photon. 

Now including higher order corrections in the energies of 
the initial photon and electrons, but still in lowest order of x%, 
instead of Iq one will obtain a more general function /. With this 
we can define an effective DC emission correction factor in the 
soft photon limit by 



(19) 



We shall now give analytic expressions for / based on the results 
obtained in the previous Section. 



3.1 .2. Direct expansion of G m 

One can also replace ^G|" v ^ in Eq. ( f20b with the more simple 
analytic expression ( fT"5l > derived in the limit of low electron tem- 
perature and energy of the initial photon. In this case the DC 
emission coefficient / can be written as the direct expansion 



(21) 



k=0 



where the definitions of the integrals /, may be found in 
Appendix iC2). The f i rst inte gral /p corre s ponds to the result 
obtained by iLightmanl dl98ll) and iThornd i 1981b in the limit 
n(v) <K 1. However, it is expected that this approximation in 
general converges very slowly and is mainly useful for simple 
analytic estimates. 

3.1.3. Alternative approach 

Above we gave the expressions obtained from the direct ex- 
pansion of I. However, this procedure results in an asymptotic 
expansion of the effective DC correction factor g™ c lt or equiva- 
lently /, which in most of the cases is expected to converge very 
slowly. Therefore, using the direct expansion ( f2Tb we again tried 
to 'guess' the correct functional form of the DC emission coef- 
ficient / and thereby to extend the applicability of this simple 
analytic expression. The general behavior of the results obtained 
in our numerical integrations showed that for e < 6 y the effec- 
tive DC correction factor g^° ft = I/Iq decreases towards higher 
electron temperatures. After several attempts we found that this 
behavior is best represented by the functional form 



^inv _ 



1 + ai 



Using Eq. (fJTJ, one finds 

21 By J x 5 ndx 
5 e J x 4 n dx 



(22a) 



(22b) 
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As will be shown below, equation d22at together with d22bt pro- 
vides a description of the low frequency DC emission coefficient, 
which for photon distributions close to a Wien spectrum is ac- 
curate to better than 5 % for temperatures up to ~ 25 keV in the 
range 0.2 < p < 50. Here we defined the ratio of the electron to 
photon temperature as p = T e /T y . For a detailed discussion we 
refer the reader to Sect. 14.21 Since this approximation involves 
only 1 -dimensional integrals over the initial photon distribution 
in numerical applications it may be more useful than the approx- 
imation ( l20b . 

3.2. Approximate inclusion of stimulated DC emission 

For broad initial photon distributions, stimulated DC emission 
arises which, in particular for situations close to full thermody- 
namic equilibrium, can become significant. Accounting for this 
effect, the DC emission integral contains additional factors of 
1 + n(vi) and 1 + n(y2) due to the presence of ambient photons 
at frequencies v\ and vi- If we are interested in the emission 
of photons at frequency V2 then 1 + n(v2) appears as a global 
factor in front of the full Boltzmann emission integral and no 
further complications arise due to this term. However, the factor 
1 + n(v\) remains inside the integrand and the dependence of the 
scattered photon frequency on the scattering angles and energies 
has to be taken into account. 

Within the Lightman-Thorne approximation, for a particular 
DC scattering event one always has vo ~ v\, since for r e = 
no Doppler boosting appears and for vo <K m e the recoil effect is 
negligible. Therefore one can write 



[1 +n(x 2 )] x 



Jr»oo 
x 4 n(x 
o 



)[1 + n(x)] dx 



(23) 



instead of Eq. ( fT8l i. In situations close to full thermodynamic 
equilibrium the main effect due to stimulated DC emission is 
only due to the factor 1 + n(xz), whereas the correction related 
to 1 + n(x) inside the integral is only of a few percent. Assuming 
n(x) = l/[e x - 1] yields I ~ 24.89 and lf m ~ 25.98/x 2 for 
X2 <k 1 which further supports this statement. 

Similarly, one can modify Eq. (l20b with the replacement 
n(x) — > n(x)[l + n(x)][l + n(xz)], but again this approach is lim- 
ited to cases when the change of the initial photon energy is very 
small (vo ~ vi) and the radiation spectrum is broad (Av/v » /3). 
This then yields 

Xoo 
x 4 n(x)[l+n(x)](G™) ih dx (24) 

Now, looking at the ratio gj™ = p tlm /[^ lm on e again expects that 
due to stimulated processes the effective DC emission correction 
factor will also be affected at the percent-level only. We con- 
firmed this statement for Planckian initial photons with T y = T e 
performing the full Boltzmann emission integral, but a more de- 
tailed treatment of this problem is beyond the scope of this paper, 
and as mentioned above below we restrict ourselves to cases in 
which stimulated DC emission is negligible. Furthermore it is 
clear that in situations when stimulated DC emission becomes 
important, DC absorption also should be taken into account and 
a full derivation of the kinetic equation for the evolution of the 
photon field under the DC process is required. In lowest order 
and for situations very close to full thermodynamic equilibrium, 
accounting for higher order corrections in the initial photon en- 
ergy and electron temperature, one can simply multiply the ki- 
netic equation for the time evolut ion of the photo n fiel d under 
DC scattering as formulated by Lightman (198l]) and iThornel 



dl98ll) . by the effective DC correction factor g s °^\ Eq. ( fT9l ), or 
alternatively g^ m = P 1 ™//^™. However, in general the situa- 
tion is much more complicated and a detail treatment of the full 
Boltzmann collision term is required, which will be left for a 
future publication. 

4. Results for Wien spectra 

4.1. Analytic expressions 

For a Wien spectrum «w(*) = No e~ x of temperature 8 y , with the 
condition A^o <k v 2 , induced terms may be neglected. Above No 
in physical units is defined by N y $ = SnNo/c^, where N y $ is the 
photon number density. Then the functions are given by 

7 ,w= 24N (25a) 

/i,w = -/o,w[21-6p]/p (25b) 

/2,w = /o,w [428.4 - 220.5 p + 15 p 2 ]/p 2 (25c) 
/a,w = -/ ,w[9141.6- 6854.4 p 

+ 1047.375 p 2 - 11.25 p 3 ]/p 3 (25d) 
h,w = / ,w[206380.8 - 205686p + 51408 p 2 

- 2480.625 p 3 - 1 1 .25 p 4 ] /p 4 . (25e) 

Here we defined the ratio of the electron to photon temperature 
as p = T e /T y . Making use of (I22al > and (I22bb one finds: 



24 N 



(26) 



w 1 + [21 - 6p] 6 y ' 
Here we replaced the electron temperature by 6 e — > p 6 y . 

4.2. Comparison with numerical results 
4.2.1. Case T e = T y 

In Fig. the results obtained for incoming Wien spectra with 
effective temperature 6 y = 6 e are shown. With increasing tem- 
perature the effective DC correction factor decreases strongly in 
both cases. This implies that for higher temperatures the effi- 
ciency of DC emission is signific antl y overes timated by the re- 
sult obtained by lLightmanl dl98ll) and lThorna (1 198 lh . For exam- 
ple, even at moderate temperatures T e ~ 4keV there is a 10% 
negative correction due to higher order corrections in the ener- 
gies of the initial photons and electrons. Since I/Io < 1, with 
Eq. ( TT3T > one can estimate that the main contribution to / has to 
come from photons with energies <±>o > y 6 e . 

In the considered case the direct formulae (f2~Tb in fourth or- 
der of the electron temperature with (l25T l are valid up to 6 e ~ 
0.05, within reasonable accuracy. As mentioned above the con- 
vergence of these asymptotic expansions is very slow. However, 
as will be shown below, the inverse formula (122ab and the cor- 
responding coefficients a for p ~ 1 provide an approximation 
which is better than 5 % up to 6 e ~ 0.2. Comparing the result 
for an incoming W ien spectrum with the approximation given 
by lSvenssonldl984l) . 



h(e e ) 



1 + 13.91 6 e + 11.05 9l + 19.92 3 



(27) 



shows that for p = 1 Eq. d27l > provides a very good fit to the 
full numerical results. The approximations based on formulae 
d20l and d26l > show a similarly good performance. However, in 
the more general case, i.e. p + 1, the approximation given by 
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Figure 5. DC soft photon production rate relative to the DC pro- 
duction rate Iq as given by the Lightman-Thorne approximation 
as a function of the electron temperature 6 e = T e /m e . All the 
curves were computed for initial photons with a Wien spectrum 
at temperature T y = T e . The full numerical results are shown 
in comparison with the direct expansion as given by equations 
( f25l ) including temperature corrections up to different orders, as 
in dicated. In additi on, the approximation formula d27T i as given 
bv lSvenssonl ([1984) is shown. The approximations based on for- 
mulae d20b and d26b show a similarly good performance but for 
clarity are not presented here. For further discussion about the 
performance of the various approximations see Sect. 14. 31 



ISvenssonl(ll984l) does not reproduce the numerical results, since 
the strong dependence of the higher order corrections on the ra- 
tio of the electron to the photon temperature was not taken into 
account. 



4.2.2. Case T e + T y 

Figure [6] illustrates the dependence of the effective DC correc- 
tion factor on the ratio of the electron to the photon temperature 
P = T e /T y for an initial Wien spectra. In general the charac- 
teristics of the ratio I/Iq can be summarized as follows: If we 
define the critical ratio p c = 7/2 ~ 2.33 of electron to photon 
temperatures, for which the first order correction to / vanishes 
(cf. Eq. ( 125bl )). then two regimes can be distinguished: (i) for 
p < p c the ratio I/Iq monotonically decreases with increasing 
electron temperature, whereas (ii) for p > p c it first increases 
and then turns into a decrease at high temperature. This suggests 
that p c separates the regions where Doppler boosting is com- 
pensating the suppression of DC emission due to higher photon 
energy. Estimating the mean photon energy weighted by Iq, i.e. 
(wo) = Q 7 J x 5 n dx/ J x 4 n dx » 5 8 y , and using ll>q ~ y 8 e for 
the condition G m (a>o,6 e ) = 1 (see Eq. (15[ ) yields p c * 7/2, 
which further supports this conclusion. 

In general one finds (see Figs. [6j> that for p < p c the inverse 
approximation ( 122at has a better performance than the direct ex- 
pansion formula (f2Tb and vice versa for the case p > p c . The di- 
rect formula performs very well for p > p c even in lower orders 
of the temperature corrections. Combining these two approxi- 
mations an accuracy of better than ~ 5% can be achieved over 
a very broad range of temperatures. If the photons and electrons 




Figure 6. DC soft photon production rate relative to the DC pro- 
duction rate Io as given by the Lightman-Thorne approximation 
as a function of the electron temperature 9 e = T e /m e . All the 
curves were computed for initial photons with a Wien spectrum 
at temperature T y and different values of p = T e /T y . The full 
numerical result, the fourth order direct formula as given by 
Eq. (|251 ) and the inverse formula from Eq. d2oT ) are given. 




Figure 7. DC emission correction factor for Wien spectra. The 
lines indicate the contours of I/Iq = const as obtained by full 
numerical integrations. Within the shaded regions the approxi- 
mation d2"0t is accurate to better than 1% and 5%, as indicated. 



have similar temperatures, then the inverse formula ( I22at is valid 
even up to kT ~ 100 keV. In Section l4~3l we discuss the range of 
applicability for the various approximations given above in more 
detail. 

In Fig|7]the dependence of the DC emission correction fac- 
tor on the temperature of the initial photons and the temperature 
of the electrons is illustrated for Wien spectra. The Lightman- 
Thorne approximation is accurate in a very limited range of 
photon energies and electron temperatures. On the other hand, 
the approximation ( |20T i works in practically the full considered 
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direct expansion formula 



E = 5% 
E = 1% 




Figure 8. Range of applicability for Wien spectra and different 
p = T e /T y : the electron temperature 6> e is shown above which 
the relative difference between the considered analytic predic- 
tion and the numerical results becomes greater than e. The nu- 
merical results were compared with the approximation (l20l >. the 
direct expansion ( |2"TT i up to fourth order in temperature and the 
inverse formula as given by Eq. 



range with very high accuracy. In contrast to the case of monoen- 
ergetic initial photons even at low temperatures of the photons 
the electrons have to be several times hotter in order to reach 
1 1 la — 1 . As discussed above this is due to the fact that the main 
emission is coming from photons with ojq > few X y . 



4.3. Range of applicability for the different approximations 

In the previous Sections we have derived different kinds of an- 
alytic approximations for the DC emission coefficient. Each has 
advantages and disadvantages. For example the approximation 
( l20t involves a 2-dimensional integral over the photon spectrum 
and the electron distribution, which is numerically more expen- 
sive than the 1 -dimensional integrals for the direct expansion 
d2TT i or the inverse approximation d22l >. One obvious advantage 
of the inverse formula d22l over the direct expansion (l2TT i is that 
only the first order corrections are needed to obtain an excel- 
lent approximation up to fairly high temperatures. Especially 
for numerical applications this is important, since higher order 
derivatives of the spectrum may lead to significant difficulties. 
However, it is necessary to investigate the range of applicability 
of each of these approximations in more detail. 

We have studied the accuracy of these approximations for 
Wien spectra more extensively. For this purpose we have com- 
pared the different approximations with the results of the full 
numerical computation and, for a given p = T e /T y , determined 
the lowest values of the electron temperature 6> e at which the ap- 
proximation deviates by more than e from the numerical results. 
In Fig. [8] we illustrate the outcome of these calculations for the 
approximation (|20T > based on the results for the monochromatic 
case, the direct expansion (l2TT i up to fourth order in temperature 
and the inverse formula from equation d22l >. The peaks in the 
curves for the direct expansion (l2lT > and approximation ( f20b ap- 
pear where the relative difference to the results of the full compu- 
tation in each case changes sign. For p < 1 the direct expansion 
breaks down very fast, while for p > 1 it works very well. The 
inverse formula works best for cases close to p — 1 . Whenever 
induced effects are negligible the approximation based on the 
results obtained for the monochromatic case is excellent, in par- 
ticular for p > 1 . 



5. Implications for the thermalization of spectral 
distortions of the CMB in the early Universe 

In the high redshift Universe (z > 2.9 x 10 5 ) injec 
tion of energy into the medium le ads to a residual u 



type distortion of the CMB tod a y dSunvaev & Zeldovic 
ll970Ullfarionov & Sunvaevl 



1975a b; Danese & de Zotti, 1982 



iBurigana. Danese & de ZottiL 11991b iHu & SilkL 119931) . In this 
case the photon occupation number can be expressed as n(v) = 
1 iy e hvikT y +n _ ^ w i];h (he dimensionless chemical potential p, 
which in general is a function of frequency v. At low frequencies 
(hv < few x lQ^kTy) the production of soft photons by the DC 
process returns the photon distribution to a Planck spectrum af- 
ter a very short time, whereas at high frequencies (hv > 0.1kT y ) 
very few photons are newly emitted and Compton scattering is 
only able to establish a Bose-Einstein spectrum with constant 
chemical potential. Therefore the time evolution of the high fre- 
quency spectrum depends critically on the rate at which soft pho- 
tons are created by the DC process and then up-scattered by elec- 
trons via the Compton process. 

The current limits on deviations of the CMB spectrum from 
a pure blackbody as obtained with Cobe ( Fixsen et all 1996; 
iFixsen & Matherl [2002) constrain the chemical potential today 
to obey \p\ < 9 x 10~ 5 , but this small residual distortion could 
have arisen from a huge energy injection at sufficiently high red- 
shifts (z > 10 6 - 10 7 ), where until today the resulting spectral 
distortions were washed out by the combined action of DC soft 
photon production and Comptonization. 
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Due to Compton scattering at these early stages the tempera- 
ture of the electrons is always very close to the Compton equilib- 
rium temperature in the given radiation field. For a /i-type distor- 
tion with id <k 1 the difference between the radiation and electron 
temperature is small and hence for estimates of the DC emission 
coefficient one may assume T e ~ T y . For small values of fi it 
is also sufficient to assume that the initial photon distribution is 
very close to a pure blackbody. As explained in Sect. 13.21 stim- 
ulated DC emission will only affect the effective DC correction 
factor at the percent level. With n(x) = 1 /[e x - 1] and Eqn. (l22l 
£d° ft ~ + 14.6 9y], which implies that at high redshifts the 
thermalization of CMB spectral distortions will be less efficient, 
since the DC photon production rate is reduced compared to the 
Lightman-Thorne formula. For a temperature of kT y ~ 4 ke V, 
higher order corrections lower the DC emissivity by ~ 10%. 
This temperature corresponds to a redshift of roughly z ~ 10 7 . 
Therefore one can expect that quantum-electrodynamical effects 
lead to significant corrections of the DC photon production rate 
for large energy injection at redshifts z > few x 10 6 . In addition, 
higher order temperature corrections to the Compton process 
may also lower the thermalization efficiency. This then will make 
the constraints placed on energy injection due to exponentially 
decaying relict particles with short lifetimes (fx < few x 10 6 s) 
tighter an d may affect the re sults obtained within the standard 
treatment (iHu & SilkL [l993a) at a level ~ 10%. However, a more 
detailed computation is beyond the scope of this work. 

Higher order corrections to the DC process are neg- 
ligible f or the theoretical comput a tions of y-type dis- 
tortions ([ Zeldovich & Sunvaev, 1969; Sunvaev & Zeld ovichL 
1970Hlllarionov & Sunvaevi|l975allbl : [Danese & de ZottiL ll982t 
Burigana. Danese & de Zottil I199U IHu & SilkL Il993l) since 
these only arise after energy release in the low redshift (z < 
10 4 ) Universe, where not only the temperature has become very 
small, but also Bremsstrahlung dominates over DC emission. 
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6. Conclusion 

Double Compton emission in an isotropic, mildly relativistic 
thermal plasma was investigated within the soft photon limit. 
Simple and accurate analytic expression for the low frequency 
DC emission coefficient have been derived, which extend the 
Lightman-Thorne approximation up to mildly relativistic tem- 
peratures of the medium and should be applicable in a broad 
range of physical situations. In particular, expressions for ini- 
tially monochromatic photons and Wien spectra were given and 
discussed in detail. 

It has been shown that the DC emissivity strongly decreases 
with higher mean energy of the initial photons leading to a sup- 
pression of the total number of newly created photons compared 
to the Lightman-Thorne approach. On the other hand increasing 
the temperature of the electrons leads to an enhancement of the 
DC emissivity. If the photons and the electrons have similar tem- 
peratures, which is the case in most physical situations close to 
full thermodynamic equilibrium, then formulae d22l) should be 
applicable up to kT ~ 100 keV with an accuracy of better than a 
few percent. Since only first order corrections are necessary for 
this approximation, it is generally most suitable for numerical 
applications. 

We also discussed possible consequences for the theoretical 
computations of the thermalization of ytz-type CMB spectral dis- 
tortions at high redshifts (z > few x 10 6 ) and found that the new 
approximations could change the results at the level of a few 
percent, possibly up to ~ 10%. However, these corrections are 
currently beyond the reach of observational possibilities. 
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Appendix A: Expansion of the relativistic Maxwell Boltzmann distribution 

The low temperature expansion of {3) leads to 



f(E) ■■ 



1 - ft. 



15 1 



345 
128 



e - + 



, (3285 345 2 15 3 25 A 
c I 1024 256 ? 16 ? 64 ? 



/ 95355 _ 3285 f 2 _ 345 f 3 _ 855 f 4 
' 1 32768 ~ 2048 256 1024 



13 f 51 f _ f f 
32 + 128 ~ 16 + 384 



(A.l) 



with f = T] /20 e and r\ = p/m s 



Appendix B: Squared matrix element for double Compton 

For the DC differential c ross section one has to calcul ate the squared matrix element \M\ 2 = e 6 X describing the DC process. This calculation was first performed by 
Mandl & Skvrme 119521) . In lJauch & Rohrlichl i 1 976 ) one can find the expression for the squared matrix element of the double Compton process (pp. 235): 



X = 2 (a b - c) [ (fl + b) (2 + x) - (ab - c) - 



4x 
~AB 



where the following abbreviations have been used 



1 1 1 
a = — H 1- — 

KQ K[ K 2 



2x\a i + b l \-2 [ab + c(l -x)]p 



(A + B) (1 +x) + jt (\-z) + 2z-(aA + bB) 2 + 



u 1 1 1 
b= — + — + — 

1 K 2 
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k a k' ki k[ ki k' 2 



(B.l) 



(B.2a) 



= KQ + Kl + K 2 



■ KQ K a + Kl K. + K2 K 2 



(B.2b) 



B = k' k[k 2 



A = K Q K\ K 2 

For the definitions of Ku k'. we used those of the original paper from|Mandl & Sk vrmel 09521) : 



-PK 
+P' ■ K 



m 2 ki = +P ■ Ki 
ml k\ = -P' ■ Ki 



A() 



+ 

K' KQ 



+ *1_ 



+PK 2 
-P' ■ Ki 



(B.2c) 



(B.3a) 
(B.3b) 



with the standard signature of the Minkowski-metric (H ). Assuming that the outgoing photon y(K 2 ) is soft compared to y(Ki ) one can expand X into orders of 

"ms of the order 0(<i> 2 2 ). Similar expan 
the results of these expansions are extremely complex and not very useful, we have omit them here. 



the frequency ai 2 and keep only the lowest order term, i.e. terms of the order 0(<i> 2 2 ). Similar expansions for the limits cjq <k 1 and B « 1 can be performed. Since 



Appendix C: Numerical solution of the Boltzmann integrals 

To solve the Boltzmann integrals we implemented two different programs, one based on the Nag routine DS1GBF, which uses an adaptive Monte-Carlo method 
to solve multidimensional integrals, the other using the Vegas routine of the Cuba LibraiyQ l Hahn, 2004). The latter turned out to be much more efficient as it 
significantly benefitted from importance sampling. The expense and performance of the calculation critically depend on the required accuracy. For most of the 
calculations presented here we chose a relative error of the order of e ~ 10~ 3 . For the integrations over different initial photon distributions we typically used 
Wmm = 10~ 4 6 y and <y max = 25 8 y . 



Integration over the electron momenta 

In order to restrict the integration region over the electron momenta for low electron temperature we determined the maximal Lorentz factor, y max , such that the 
change in the normalization of the electron distribution was less than a fraction of the required accuracy. Instead of p we used the variable £ = rj 1 /2 8 t with r\ = p/m e . 
For high electron temperature it turned out to be more efficient to use the normalization of the electron distribution itself as a variable. For this we defined 



N(x) = m { 



y 2 - i/'(y"ie)dr, 



(C.l) 



which for x — » oo becomes N = N e . In actual calculations one has to invert this equation to find the function x(N). This was done numerically before the integrations 
were performed for a sufficiently large number of points (n ~ 512) such that the function x(N) during the integrations could be accurately represented via spline 
interpolation. Here it is important to bear in mind, that N(x) is rather steep for x ~ and x ~ 1. 



4 Download of the Cuba Library available at: http://www.feynarts.de/cuba/ 
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C.1. Definition of h 

Inserting the direct expansion of G m in low energies of the initial photon and electron temperature Eq. il5i into Eq. )20t and collecting the different orders in ft c , one 
can define the functions /, of equation J2U as 



/o = £>4 

21 6y 



(C.2a) 
(C.2b) 



441 0y 357 

h = 15 £> 4 D 5 + 

10 e e 25 



4 40 6». 5 25 



2 



7618 
175 



45 



3969 



8 ft 



8568 



£>6 



34281 
35 



2>7 

3 



£>7 + 



21498 
175 



(C.2c) 
(C.2d) 
(C.2e) 



where we introduce the integrals £); as 



2) ; = J x'ndx 



(C.3a) 



over the initial photon distribution n(v). 



